Cosmological aspects of inflation in a supersymmetric axion model 



O 

(N 

C 

(N 



43 

i , 

Or 
CD 

43 



> 
(N 

(N 

o 



X 



Masahiro Kawasaki, 1,2 Naoya Kitajima, 1 and Kazunori Nakayama 3,4 

1 Institute for Cosmic Ray Research, University of Tokyo, Kashiwa, Chiba 277-8582, Japan 
Institute for the Physics and Mathematics of the Universe, 
University of Tokyo, Kashiwa, Chiba 277-8568, Japan 
1 Department of Physics, University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
^KEK Theory Center, Institute of Particle and Nuclear Studies, KEK, Tsukuba, Ibaraki 305-0801, Japan 

(Dated: January 19, 2013) 

We show that the hybrid inflation is naturally realized in the framework of a supersymmetric axion 
model, which is consistent with the WMAP observation if the Peccei-Quinn symmetry breaking scale 
is around 10 15 GeV. By solving the post-inflationary scalar dynamics, it is found that the scalar 
partner of the axion, saxion, oscillates with large amplitude and its decay produces a huge entropy 
and dilutes the axion. As a result, the axion coherent oscillation can be the dominant component 
of the dark matter in the Universe. Cosmological gravitino and axino problems are solved. 



I. INTRODUCTION 

The standard model of particle physics has succeeded 
in explaining many experimental facts until now, but 
there exist several problems which have not been solved 
yet. One of them is the strong CP problem. The La- 
grangian of quantum chromodynamics (QCD) allows the 
CP-violating term, but the neutron electric dipole mea- 
surement implies that QCD preserves CP and hence 
the CP-violating term is stringently constrained. The 
most popular solution for the problem was proposed by 
Peccei and Quinn [l|. They introduced an additional 
global U(l) symmetry, called Peccei-Quinn (PQ) sym- 
metry, written as £/(1)pq. When the PQ symmetry is 
broken spontaneously, a pscudo-Nambu-Goldstonc boson 
appears, which is called an axion. The axion acquires its 
mass through the QCD instanton effect and the axion set- 
tles down to the potential minimum (= vacuum) where 
CP is preserved @|. 

Another well-known problem of the standard model is 
the gauge hierarchy problem. The electroweak scale is 
unstable against the radiative correction. The squared 
Higgs mass receives a quadratically divergent correction 
which is thought to be, say, the grand unification scale 
~ 10 16 GcV. Thus keeping the Higgs mass around the 
weak scale requires unnatural fine-tuning. This prob- 
lem is solved if we introduce the supersymmetry (SUSY) 
which is the symmetry with respect to the replacement of 
bosons with fcrmions. In SUSY theory, the quadratically 
divergent quantum correction is canceled [3[. Thus we 
have a good motivation for considering the axion model 
in the framework of SUSY. 

On the other hand, the cosmological observations sup- 
port that the Universe underwent an accelerated expan- 
sion, called inflation, at the very early stage of the Uni- 
verse [|| . It is generally considered that inflation is driven 
by potential energy of a scalar field, inflaton. Thanks to 
inflation, problems of the standard cosmology such as the 
horizon problem and the monopole problem are solved. 
In addition, inflation can produce the density perturba- 
tion which accounts for the large scale structure of the 



present Universe. The observation by the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) satellite confirmed 
the nearly scale invariant density perturbation, which is 
a generic prediction of inflation. 

The cosmological observations also revealed the ex- 
istence of the nonbaryonic dark matter (DM), which 
occupies 23% of the total energy of the present Uni- 
verse- [|. The standard model (SM) in particle physics 
does not contain candidates for DM, and hence we need 
the physics beyond the SM in order to explain the exis- 
tence of DM. 

In the previous paper Q , we pointed out that a SUSY 
axion model naturally causes hybrid inflation and the 
axion becomes the dominant component of DM. In the 
hybrid inflation model, the PQ scalar fields play a role 
of the waterfall field and the PQ phase transition takes 
place at the end of inflation. To account for the ampli- 
tude of the density perturbation observed by the WMAP, 
the PQ scale f a must be of order of 10 15 GeV. This value 
seems to be too large because there is an upper bound 
fa ^5 10 12 GeV in order for the axion not to overdose the 
Universe. However, in the SUSY axion model it is found 
that the axion is diluted by the late-time entropy produc- 
tion via saxion decay and the upper bound is relaxed to 
10 15 GeV HQ , which is consistent with the condition for 
the correct density perturbation. In this paper we fur- 
ther describe the details of the scalar field dynamics after 
inflation in this model by performing detailed numerical 
calculations and we show that the saxion field deviates 
from the true minimum by the finitc-tcmpcraturc effect 
and the saxion oscillation has a large amplitude, which 
leads to a huge entropy production and dilutes the axion 
at a later epoch. We also study several problems in our 
model such as the axion domain wall and baryogenesis 
and propose several solutions to these problems. 

This paper is organized as follows. In Sec. 2 we briefly 
introduce a SUSY axion model. In Sec. 3 we review the 
SUSY hybrid inflation model and show the consistency 
with observation. Particularly we point out that the su- 
perpotential of the SUSY hybrid inflation model has the 
same form as that of the SUSY axion model, and the 
PQ scale is determined to be ~ 10 15 GeV for the PQ sec- 
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tor to induce a hybrid inflation and correctly reproduce 
observed density perturbation. In Sec. 4 the dynamics 
after inflation is discussed and we see that the saxion os- 
cillation with large initial amplitude is induced. In Sec. 5 
it is shown that the late-time entropy production due to 
the saxion decay necessarily takes place, and as a result, 
the axion coherent oscillation can be the dominant com- 
ponent of DM. In Sec. 6 we discuss the fate of topological 
defects, such as axionic strings and domain walls. In Sec. 
7 we present a mechanism to create a correct amount of 
baryon asymmetry under the late-time entropy produc- 
tion. In Sec. 8 a variant type of SUSY axion model is 
presented, which causes a so-called smooth-hybrid infla- 
tion and we describe some cosmological aspects of the 
model. We conclude in Sec. 9. 



II. SUPERSYMMETRIC AXION MODEL 
A. The potential of the SUSY axion model 

Let us describe the SUSY axion model. Here we as- 
sume the gravity-mediated SUSY breaking. The super- 
potential for the SUSY axion model is given by 



W = - fl) + A$II, 



(1) 



where, S is a gauge singlet superfield and has a zero 
PQ charge, and ^ and ^> are the PQ superfields that 
are gauge singlets and have +1 and — 1 PQ charges, re- 
spectively. The PQ fields contain the axion (a), saxion 
(a, the scalar partner of the axion), and axino (a, the 
fermionic supcrpartner of the axion). Here f a is the PQ 
symmetry-breaking scale and k is a dimensionless cou- 
pling constant assumed to be real and positive. X(X) is 
the superfield interacting with a PQ field at tree level and 
has some PQ charges as well as gauge charges through 
which it interacts with the minimal supcrsymmctric stan- 
dard model (MSSM) fields. The superpotcntial also has 
an i?-symmctry. The charge assignments of the fields in 
the present model are shown in Table |TJ In particular, 
for the Kim-Shifman-Vainshtein-Zakharov (KSVZ) (or 
hadronic) axion model [§| , X and X are additional heavy 
quarks, denoted by Q and Q, that have color charges. 
For the Dine-Fischlcr-Srednicki-Zhitnitsky (DFSZ) axion 
model @, X and X are identified as MSSM Higgscs, H u 
Id- 



TABLE I: Charge assignments on the field content 
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According to the superpotcntial (JXJ) , the F-term scalar 
potential is derived as 



- / q t + n 2 \s\ 2 m 2 + (2) 



where we denote the scalar fields using the same symbol 
as the superfields and we set X and X to be zero as- 
suming that they have large Hubble masses and quickly 
settle down to zero during inflation. The global minimum 
of this potential is located at S = and tyfy = f%. Here 
it should be noted that there exists a flat direction along 
which the scalar fields do not feel the potential, ensured 
by the U(1)pq symmetry extended to a complex U(l) 
due to the holomorphy of the superpotcntial [10| . The 
flat direction is lifted up by the SUSY-breaking effect, 
leading to the following soft SUSY-breaking mass terms 



U so ft = ciTOg /2 |*| 2 + c 2 m\ /2 



(3) 



where m.3/2 is the gravitino mass and c\ and c 2 are 
real-valued constants that are positive and order unities. 
With this soft SUSY-breaking potential, the radial com- 
ponents of the PQ fields \^>\ and |^| are stabilized at 
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respectively. 2 The saxion field a is defined by the devia- 
tion of l^l from the vacuum expectation value ^ along 
the flat direction. 

Near the vacuum expectation values ((3]), the axion a 
and saxion a are related to the PQ fields as 



v exp 



x a 



V2F a 



v exp 



V2F a 



(5) 



where F a is determined by requiring that a and a are 



canonically normalized and given by F a = \J v 2 



1 In the DFSZ model, the coupling constant A must be very small, 
say, A ~ 10~ 12 for f a ~ 10 15 GeV, in order to produce a sizable 
/^-term. This might be a tuning, but it is relaxed by changing 
the relative PQ charge assignments between and H u (Hd). 

For example, if the PQ charges of H u and Hd are — n, where n(> 
1) is a positive integer, the allowed term in the superpotcntial 
is \<b 2n H u H d /M 2n ~ 1 with some cutoff scale M, which might 
be the Planck scale. In this case the amount of tuning for the 
coupling constant A is relaxed. The phenomenology discussed in 
the following sections is not modified by the choice of PQ charges 
for H u and Hj. 



B. The decay of the saxion 

In this subsection, we derive the decay rate of the sax- 
ion which is important in the later section. The kinetic 



2 See Refs. [l l!4l7t ] for other types of the saxion stabilization mech- 
anisms and their cosmological issues. 
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terms of the PQ fields lead the interaction of the saxion 
with the axion [3, as 



IT 



(6) 

where £ is defined as £ = (w 2 — v 2 )/F 2 , which is generally 
of order unity unless v ~ w i.e. ci ~ ca. From @, the 
decay rate of the saxion into two axions is estimated as 



superpotential and Kahlcr potential which are responsi- 
ble for inflation are written as 



K 



m 2 



\s\ 2 - 



fci 



Ml 



M 2 



kss-^ttt + 



6MJ, 



IU = K S(**-/ a 2 ) + Wo, 



(10) 
(11) 



64tt F 2 



(7) 



Note that the decay rate into two axions can be sup- 
pressed by tuning C\ and C2 as c\ ~ C2- 

In the KSVZ axion model, the dominant decay mode 
into the SM particles is that into two gluons. The decay 
rate is calculated as 



>99 



32tt 3 F 2 



(8) 



where m a is the mass of the saxion and a s is the QCD 
gauge coupling constant. Comparing flSJ with (|7J|, the 
dominant decay mode of the saxion is generally the two- 
axion decay in the KSVZ axion model unless c\ ~ C2. 
In a minimal setup where the SUSY-brcaking masses for 
^/ and only come from the supergravity effects, c\ is 
equal to C2 at tree level. The radiative correction through 
heavy quarks, which is relevant only for vff, produces a 
difference between c\ and C2 at the low energy scale [l8[ . 
However, even if such effects are taken into account, the 
two-gluon decay can be the dominant mode. 

In the DFSZ axion model, there exists a tree level cou- 
pling of the saxion with the standard model Higgses, as 
£ int = A 2 cr 2 (|#°| 2 + \H a d \ 2 ). Thus the saxion decays into 
Higgses with decay rate given by 



1 ra a 
8^7f 



4m 2 



1/2 



(9) 



where S and ^(^/) take roles of the inflaton and water- 
fall fields, respectively; Wo = m^i^Ml is the constant 
term which cancels the vacuum energy coming from the 
SUSY-brcaking effect and makes the cosmological con- 
stant nearly zero in the present Universe; and ks, ki, k 2 , 
and kss are dimensionless real-valued coefficients. As 
seen later, the PQ scale f a sets the energy scale of infla- 
tion. 

In the framework of the supergravity, the scalar poten- 
tial is obtained from 



V 



,K/M 2 f 



K ij DA¥D*W* - 3 



\W\< 



M 2 P 



(12) 



where Mp is the Planck scale, DiW = W{ 



KiW/Mf,, 

Kij*, and subscript i means the derivative with 
respect to the field <fo = {S, In the global SUSY 

limit, i.e. Mp — > oo, the F-term potential that comes 
from the superpotential (fTTj) becomes 



V = K 2 \W - f 2 \ 2 + K 2 |S| 2 (|#| 2 + |^| 2 ). 



(13) 



The minimum of the potential lies at S — and "I"!' = J 2 , 
where V = 0. However, if \S\ takes a large value, the 
(local) minimum appears at ^ = ^ = where V = K 2 f^. 
Since the potential for |5| > f a with * = * = is flat, 
inflation takes place there. Once |5| rolls down from a 
larger value and reaches |5| = / Q , the waterfall behavior 
turns on and the inflation ends suddenly. To see this, we 
rewrite the mass terms for \t and \& in the potential (JT3J) 
as 



where \i = A(W) gives the Higgsino mass. In the DFSZ 
axion model, this decay mode may dominate over the 
two-axion decay. 



III. HYBRID INFLATION IN A SUSY AXION 
MODEL 



In this section, we show that the axion model described 
in the previous section is also considered as a SUSY hy- 
brid inflation model (2(| [2l|. 3 The relevant parts of the 



V 




(14) 



So there are two mass eigenstate with squared masses 
k 2 (|S'| 2 ± / 2 ). Therefore when \S\ < f a , one of the mass 
cigenstates becomes tachyonic and gets a nonzero vac- 
uum expectation value. After that, S approaches to zero 
and $ and \& relax to the flat direction given by *ff$! = J 2 . 

We rewrite the inflaton S as S = ipe 16 /v2. During the 
inflation period, including supergravity effect, the scalar 
potential is given by 



3 The different SUSY hybrid inflation models related to PQ sym- 
metry are found in 122 . [2311 . 



V = K?ft 1-fc 



'2M 2 P 



is 



9 



8M P 



(15) 



+2 K f 2 m 3/2 (S + S*)+ AVi_ loop , 



4 



where 75 = 1 — Tks/2 + 2k 2 s — 3kss an d AVi_i oop is 
the one-loop radiative correction. Note that the squared 
masses of the scalars and fermions in superfields S, "J, 
and are 

S # 

scalar : x 2, n 2 (\S\ 2 ± M 2 ) x 2 (16) 

fermion : x 2, k 2 \S\ 2 x 4. (17) 

Then, the one-loop radiative correction in the scalar po- 
tential [HI is calculated as 

A^ 1o o P = E(-1) F 7^1^ = #^)> (18) 

i 

where the sum is taken over the field degrees of freedom 



and F = for scalar and F = 1 for fermion. Here A is 
some cutoff scale, x is defined as x = \S\/ f a = (p/y/2f a , 
and 



F(x) 



(.t 4 + 1) In 



x 



1 



2x 2 ln 



c 2 + 1 



1 



+ 21n- 



k 2 M 2 x 2 



A 2 



(19) 



In particular, for x 3> 1, i.e. (/? ^> \[2f a . the one- loop 
correction can be approximated as 



AVi_ 



loop 



16tt 2 



2A 2 



(20) 



From the scalar potential, the slow-roll parameters [25| 
are derived as 



J 



Mp 



^ = -k. 



V2M F 



is 



2V2MI 



1 Mf 



167T 2 / a 



F'{x) + 



2m^/2Mp cos 9 



(21) 



V" 

7 1 =M 2 P — = -k s 



37£ 



<P 



2M 2 p 



^ + F» {X ), 



167T 2 / 2 



(22) 



where the prime denotes the derivative with respect to 
ip. The number of e-foldings is calculated from 

N = [, Hdi = -k[M dv ' (23) 

where the subscript / means the moment when the in- 
flation ends. The number of e-foldings after a comoving 
scale k leaves the horizon is written as 



N e =56 + In 



0.002 Mpc" 1 



k 



i) 



-1 f Tr \ + h ( Hl 

3 n I 10 10 GcV J 3 n I 10 13 GeV 



(24) 



where Tr is the reheating temperature and Hi is the 
Hubble constant during inflation. Using these quantities, 
the power spectrum of the curvature perturbation and 
the scalar spectral index n s , can be calculated as follows: 



Pt 



1 



V"' 



l2ir 2 M P V' 2 ' 



n s ~ 1 - 6e + 2?/. (25) 



We can compare these quantities with the results of the 
WMAP observation d, 



P C = (2.441±° ^ 8 2 ) x 10- 



ns = 0.963 ±0.012. (26) 



According to [2l| , to explain the correct magnitude of 
the density perturbation, f a ~ 10 15 GeV is needed and 
by adjusting the coupling constant k and the coefficient 
of a nonminimal Kahler term k$ as k ~ 10" 3 and k s ~ 
10 -2 , the scalar spectral index fits well with the WMAP 
observation. 



IV. THE POST-INFLATIONARY DYNAMICS 

In the previous section, we show that the successful 
inflation occurs in the SUSY axion model if the PQ scale 
is around 10 15 GeV. This value seems to be too large since 
fa ^$ 10 12 GeV is needed in order for the axion not to 
overdose the Universe. Note that it is not allowed to tune 
the amplitude of the axion in the present model, since 
the PQ symmetry is restored during inflation and the 
misalignment angle takes random values for each small 
patch of the Universe after inflation. However, if the 
late-time entropy production takes place, the axion can 
be diluted so that the bound can be evaded 0- In this 
and the next section, we discuss the dynamics of the PQ 
scalar fields after inflation and show that the saxion starts 
to oscillate with large initial amplitude and it results in 
a huge entropy production through its decay. 

Soon after inflation ends, the inflaton S begins to oscil- 
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isi/f a - 
m/f a 







10 3 10 4 10 5 

f a x time 



FIG. 1: Time evolution of S (thin red line), ^ (thick green 
line), and ^ (thick blue line) after inflation, as a function of 
time. We have taken f a = 0.1A1p,k = 0.01. (The green line 
and the blue one coincide.) 




f, x time 



FIG. 2: Evolution of the saxion a, inflaton S, and Hubble 
parameter H. We have taken f a — 0.1A/p,« = 0.01 and 
Cl = 1, C2 = 4. 



late around its vacuum expectation value (VEV), S ~ 0. 4 
Since the inflaton oscillation behaves as matter, its ampli- 
tude decreases inversely proportional to the cosmic time, 
i.e. oc 1/t oc a~ 3 / 2 . The waterfall fields quickly roll down 
to the flat direction, vE"^ = f%, as shown in FigJT] In 
fact, the flat direction is not flat at this stage, because 
both PQ fields obtain masses of k\S\ before S decays. 
Thus, the PQ fields are stabilized at = ^> = f a at this 
stage and both fields oscillate around f a . 

The evolution of the inflaton and the saxion is shown 
in Fig. [5] After inflation, PQ scalars oscillate around 



4 Precisely speaking, the VEV of S deviates from zero because of 
the supergravity effect. The second term in (15) makes (S) ~ 
m 3/2/ K a t the true minimum. 



|^| = 1 4' | = f a because they have masses of k\S\. When 
the gravitino mass exceeds k\S\, the saxion begins to 
move toward the true minimum ([4]). However, the Hub- 
ble parameter is already smaller than the gravitino mass 
and hence the friction is not efficient. As a result, the sax- 
ion adiabatically approaches to the true minimum with- 
out oscillation. A key point is the existence of the mass 
of k\S\ for PQ scalars, which is larger than the Hubble 
scale. Thus the saxion oscillation is not induced [2(|, and 
hence no late-time entropy production takes place unless 
S decays before the Hubble parameter becomes equal to 

The inflaton S can also decay into the axino pair, which 
is represented by the combination (ip — i/j)/-\/2, where 
ip and ip denote the fermionic component of *S> and \&, 
respectively. This process, however, cannot be used as 
a reheating process because the produced axinos are not 
thermalized. 

Therefore, the inflaton is required to decay soon after 
inflation in the present scenario. The inflaton can decay 
through the following superpotential: 

W = kSYY, (27) 

where Y(Y) is the chiral superfield which has some gauge 
charges but no PQ charges and A: is a coupling constant 
assumed to be real and positive. For example, in the 
framework of the KSVZ axion model, it can be the MSSM 
Higgses, Y — H u and Y = Hd- 5 This term also generates 
the sizable /i-tcrm for k ~ k since, as noted earlier, the 
VEV of S is given by (S) ~ m 3 / 2 /«;. Then, the decay rate 
is given by ~ l/(87r)fc 2 TOs and the reheating temperature 
after inflation, Tr, is estimated as 

(28) 

taking mj = nf a into account. 6 Thus the reheating tem- 
perature can be about 10 11 GeV for k = k = 10~ 3 and 
fa = 10 15 GeV. This is high enough to induce the saxion 
oscillation, as seen in the following. 

One complexity arises from thermal effects on the 
scalar potential. The reheating process thermalizes the 
Universe and the decay products of the inflaton inter- 
acts among others in the thermal bath. In particular, in 
the KSVZ model, the * interacts with heavy quarks Q 
and Q, and they interact with MSSM particles through 



5 The inflaton can also decay into gluons through the heavy quark 
loops in the KSVZ model, but it is subdominant under the pres- 
ence of the term ]27t . 

6 The introduction of the superpotential II27I I may modify the 
inflaton dynamics significantly. First, if k > K, the additional 
Coleman- Weinberg correction larger than the original one arises. 
Second, if k < K, the inflation ends at the point where Y(Y) 
becomes tachyonic. Here we assume k ~ k so that the original 
inflationary dynamics is not much affected. 
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f a x time 



FIG. 3: Evolution of the PQ fields, the Hubble parameter and 
the thermal mass for « = 0.01, /„ = O.OlAfp, and T R = 0.1/ o . 
The soft SUSY-breaking mass is not included. 
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FIG. 4: Evolution of the PQ fields, the Hubble parameter, and 
the thermal mass for m 3/2 = KT 6 /a, k = 0.01, /„ = O.OIMp, 
and Tr = O.l/o. 



the QCD couplings. 7 Then the two-loop effect induces 
a finite-temperature correction on the potential of \1/ ex- 
pressed as |27| 



V t 



th 



a?T 4 In 



1*1! 

T 2 



(29) 



This thermal-log potential lifts up the flat direction, and 
|*| (| 'J |) tends to roll down to a smaller (larger) value. 
With use of the condition for the flat direction 'F'F = f 2 , 
1^1 obtains an effective thermal mass from the thermal- 
log potential ([29]) 



a.s 



'th 



ajT 4 

W 2 



and the equation of motion of * is written as 
* + 3H^> + m 2 h ® = 0. 



(30) 



(31) 



The evolution of * is shown in Fig. [3] Just after the re- 
heating completes, both and ^ sits at f a as explained 
earlier, and hence the thermal mass is larger than the 
Hubble parameter at this stage. Thus <F rolls down the 
thermal potential and \^>\ increases until the effective 
mass (|3"0")) becomes comparable to the Hubble parame- 
ter. Then the "J stops rolling and gets frozen. The frozen 
value of I M/ 1 is determined at the time H ~ mth and es- 
timated as 



1*1 



(32) 



After the temperature decreases and soft SUSY- 
breaking masses (~ TO3/2) exceeds the thermal mass, the 



true minimum (0]) appears. On the other hand the PQ 
fields still remain frozen at (|32|) until the Hubble pa- 
rameter becomes smaller than the soft SUSY-breaking 
masses. Finally, after that, the PQ fields restart oscilla- 
tion around its vacuum expectation value along the fiat 
direction. After all, the saxion begins to oscillate with 
an initial amplitude 



Oi ~ a s Mp, 



(33) 



at H ~ 7713/2- These features are confirmed by numerical 
calculation and seen in Fig. SJ Although we have chosen 
specific parameters for a numerical reason in the figure, 
qualitative arguments do not change for more realistic 
parameter choices. We will see in the next section that 
the decay of the saxion dilutes the axion and its present 
cosmic density is acceptable. 



V. THE LATE-TIME ENTROPY PRODUCTION 

In this section, the cosmological scenario after the sax- 
ion begins to oscillate is described. As previously noted, 
due to the finite-temperature effect after reheating, the 
initial amplitude of the saxion oscillation is as large as 
CTj ~ a s Mp. The saxion begins to oscillate when the 
Hubble parameter becomes equal to the saxion mass, m a , 
in the radiation dominated Universe. At that time, the 
ratio of the energy density of the saxion to the entropy 
density, which is a constant until the saxion decay, is 
calculated as 



7 In the DFSZ model, we also have to introduce such heavy quarks 
in order to solve the domain wall problem. This will be discussed 
in Sec. 6. Therefore, the presence of the finite-temperature po- 
tential is generic for both the KSVZ and DFSZ models. 
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-T, 



M P 



8 

4 x 10 6 GeV 



VlTcV/ 



1/2 



a 8 Mp 



(34) 
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where the subscript i represents the value at the begin- 
ning of the saxion oscillation and 



90 



1/4 



y/M. 



pm a ~ 2 x 10 1U GeV 



VlTeV/ 



1/2 



(35) 

Here we have used the value <?* = 228.75 as the rela- 
tivistic degrees of freedom. Since the saxion coherent 
oscillation behaves as a matter, the energy density of the 
saxion soon dominates the Universe if the initial ampli- 
tude is sufficiently large. 

Here we estimate the saxion decay temperature. From 
the saxion decay rate given previously ©, the saxion 
decay temperature is estimated as 



T a ~ 5 MeV 



x 3/2 

10 TeV/ 



10 15 GeV 

fa 



(36) 



for the KSVZ axion model. In order not to destroy the 
success of the Big Bang nucleosynthesis (BBN), the con- 
dition To- > a few MeV must be imposed 30] . Here we 
have assumed that the decay into axions are suppressed. 
As noted earlier, this is actually the case if c\ ~ C2. 8 
Note also that rather large saxion mass of ~ (lOTeV) 
is required for successful reheating. If the MSSM parti- 
cle masses are around ITeV, the saxion also decays into 
them with a similar decay rate to the SM particles. This 
is problematic since the lightest SUSY particles (LSP) are 
nonthermally produced and its abundance is too large to 
be consistent with the WMAP observation, if the LSP 
is stable. Thus we need to either introduce an i?-parity 
violation in order for the LSP to decay rapidly [31(, or 
assume all SUSY particles are so heavy that the saxion 
cannot decay into them kincmatically. 

In the case of DFSZ axion model, from the decay 
rate J9]), we obtain 



T a ~ 5 MeV 



m a \V2 /10 15 GeV\ / fi x 



ITeV/ V fa 



(37) 



In this case T a > 1 MeV is satisfied for /„ ~ 10 15 GeV 
and m a ~ 1 TeV. Therefore, the saxion decay into SUSY 
particles can naturally be forbidden for TeV scale SUSY 
masses. 

We can then calculate the dilution factor defined as the 
ratio of the entropy density before and after the saxion 
decay as 

_ gbefore _ 3 / p a (t j) ^ 1 
Saftcr 4 \ 



2 x 10 



, T a WlTcV 
lMeVy 



1/2 / „ r \ 2 

' a s Mp 



(38) 



The nonthermal axions produced by the saxion decay con- 
tributes to the extra radiation of the universe. If the branch- 
ing ratio into the axions is around 0.1, it may account for the 
recently claimed existence of the extra radiation 0, l28t |29|| . 



A. The axion abundance 

Now let us estimate the axion abundance under the 
entropy production after the saxion decay. The axion 
starts to oscillate when H(T) = m a (T), where m a (T) is 
the temperature-dependent axion mass given by [32|, |33[ 

m (T) ~ {°- 08TO a( A QCD/T) 3 - 7 for T > Aqcd/tt 
\m a for T < Aqcd/tt, 

(39) 

and m a is the zero temperature mass of the axion given 
by 



m a ~ 6 x 10 _9 eV 



10 15 GeV 

fa 



(40) 



The energy density of the axion is given by p a = 
(l/2)rn^f^9l, where 6\ is the initial misalignment an- 
gle. Here and hereafter the subscript "1" refers to the 
time when the axion starts to oscillate. Because the PQ 
symmetry is broken at the end of the inflation, the initial 
misalignment angle takes random values in each horizon. 
Therefore, we take the averaged initial misalignment an- 
gle as 7T /y/3. The ratio of energy density of the axion to 
that of the saxion at the beginning of the axion oscillation 
is given by 



Pa(*l) _ * 2 fl 
Ar(*l) 



18M p ■ 



Then the axion density to entropy ratio is 

Pa(h) 7T 2 / Q 2 Pa(h) 



Si 



18Mp 



si 



(41) 



(42) 



Using Eq. (|38[) . the present axion density to entropy ratio 
is written as 



Pa(*0) _ 7r l (T ya(«l 



•ST) 



si 2AM 2 P ^{T l y 



(43) 



where ^ is defined as the ratio of the temperature- 
dependent mass to zero temperature mass of the axion 
: £(Ti) = m a (Ti)/m a . For the PQ scale /„ as large 
as 10 15 GeV, the finite-temperature effect is negligible at 
t = t\ and hence it is approximated as £(Ti) « 1 Q. 
To show this, using the scale factor dependence of the 
energy density of radiation p r cx a~ 3 / 2 when the saxion 
decays into radiation gradually in the saxion-dominated 
Universe, the relation 



P<r(t<j) 



(44) 



is derived. Using p a (t a ) = iT 2 a M pl p a {h) = 
3m 2 (Ti)M|,, we get 



m a (Ti) 



1/4 



(45) 
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and hence the temperature that the axion begins to os- 
cillate is estimated as 



axino abundance from thermal production [40[ , after the 
dilution, is estimated as 



r 1 c0.2GeV( 1Ql ^ eV ) 



0.13 



T 

-*- IT 



1 McV 



(46) 



Thus, for f a ~ 10 15 GeV, the axion starts to oscillate be- 
low the QCD scale and then the finite-temperature effect 
can be negligible. 

Finally, making use of the present ratio of the crit- 
ical density to the entropy density /9 C r,o/so — 3.64 x 
10 _9 /i 2 GeV, the density parameter of the axion is calcu- 
lated as 

Therefore, the axion can be the dominant component of 
the DM by taking account of the entropy production pro- 
cess from the saxion decay. As described later, the ax- 
ionic strings emit the axion and its contribution is com- 
parable to the coherent oscillation one. 



B. The gravitino and axino abundance 

It is well-known that in a SUSY theory, gravitinos are 
efficiently produced at the reheating and they may be 
cosmologically harmful depending on the reheating tem- 
perature. In the present model, however, the late-time 
entropy production dilutes the harmful relics. Let us see 
it. 

Gravitinos are produced both thermally [34j and non- 
thermally [35[ from the inflaton decay, but diluted suffi- 
ciently. The thermally produced gravitino abundance, in 
terms of the number to entropy ratio, is estimated as 



c4 x io-- 



10 n GcV/ VlMcV 



lTeVA 1/2 (a s M P 



(48) 



The nonthermally produced one through the inflaton de- 
cay is given by 



T R J \lMeV 



1/2 



fa 



10 15 GcV/ U0 



k \ 2 ( a s Mp 



m 

2 



(49) 



These satisfy the bound on the unstable gravitino abun- 
dance from BBN m 3/2 Y 3/2 < lCT 13 -l(r 9 GeV for m 3/2 ~ 

i-io Tev mm. 

The axino, which is the fermionic supcrpartner of the 
axion, might also have significant effects on cosmology 
once they are produced in the early universe [38U4l| . The 



K (TP) ~1 x 10 



,TTeV\ 1/2 / T R \ / T a \ 

n a ) vio n Gevy viMcvy 



10 15 GeV\ fa s M P 



fa J V G 



(50) 



In the present model, the axino mass comes from the 
VEV of S. It generates the axino mass of m a = k(S) ~ 
TO 3 / 2 - Thus the axino mass is comparable to the grav- 
itino. If the axino is not the LSP, it has a similar lifetime 
to the saxion in the KSVZ model, and it decays before 
BBN. The constraint is given as Y a < 10~ 12 so as not 
to produce too much LSPs. If the axino is the LSP, the 
bound reads m a Y a <4x 10~ 10 GcV. In both cases, the 
constraint is satisfied as is seen in Eq. (|51[) . 

A more important contribution to the axino comes 
from the direct decay of the inflaton S. As noted in 
Sec. IIV1 the S decays into the axino pair. The branch- 
ing ratio into the axino pair, B a , is evaluated as B a ~ 
K 2 /(4fc 2 ). The axino abundance produced in this way is 
written as 



K (NTP) ~1 x 10- 12 i?a 



ITcV 



1/2 



T a \ / !0 15 GeV \ 

lMcvJ { fa J 



10 n GeV, 
10" 3 \ fa,M F 



(51) 



This is close to the observational upper limits shown 
above. 9 

Therefore, the cosmological problems related to the 
gravitino and axino are solved due to the late-time en- 
tropy production from the saxion decay. 



VI. AXIONIC STRINGS AND DOMAIN WALLS 

After the PQ phase transition, the /7(1)pq symmetry 
is spontaneously broken, and axionic strings arc formed 
accordingly through the Kibble mechanism. The explicit 
breaking of the C/(1)pq by the anomaly effect further 
breaks it down to a Zn symmetry, where N is the color 
anomaly number. At the QCD phase transition, Zn sym- 
metry is spontaneously broken and domain walls are pro- 
duced. Therefore, a complicated network of topological 
defects appears such that each string is attached by N 
domain walls. In particular, domain walls are harmful 
because they soon come to dominate the Universe. 



9 The branching ratio Bg can be suppressed if k > k is allowed. 
However, the condition k > k affects the inflationary dynamics 
as shown in footnote [6] In this case, the larger f a is derived from 
the CMB normalization, which is not desirable. 
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In the case of N = 1, however, there is no such a 
domain wall problem. The reason is following. At for- 
mation, the energy of such a domain wall-string system 
is dominated by the string but soon the energy of the do- 
main wall becomes dominant. Once domain walls dom- 
inate the wall-string system, the surface area tends to 
become smaller by the surface tension of the domain 
walls. As a result, the domain wall shrinks and disap- 
pears. Thus the cosmological problem does not arise (42J. 
Through these processes, axionic strings and walls lose 
their energy emitting axions. The emitted axion obtains 
a mass after the QCD phase transition and gives signif- 
icant contributions to the present energy density of the 
Universe in addition to the coherent oscillation. Accord- 
ing to Refs. [H, HH, the energy density of the axion ra- 
diated from strings is comparable to that of the coherent 
oscillation. This is true even in the presence of late-time 
entropy production, since both the coherent oscillation 
and string-induced axion are diluted in a similar way. 
In the KSVZ model with one set of heavy quarks, the 
color anomaly number is equal to 1 and no domain wall 
problem arises. 

In the case of N > 2, however, the produced domain 
walls have complicated structure. N sheets of walls are 
attached to each string. This wall-string system domi- 
nates the Universe immediately and the present energy 
density of the domain walls becomes much larger than 
the critical energy density of the Universe. It is a cos- 
mological disaster. This is the case for the DFSZ axion 
model. One of the solutions to the problem is to intro- 
duce some heavy quarks like the KSVZ axion model in 
order to make the color anomaly number 1 [45|. Actu- 
ally these heavy quarks, introduced to solve the domain 
wall problem, serve as a source for the finite-temperature 
potential which induces the saxion oscillation with large 
amplitude, as discussed earlier. 



VII. BARYOGENESIS THROUGH THE 
AFFLECK-DINE MECHANISM 



where M is some cutoff scale. Then the potential for the 
AD field is written as 



Vs(<t>) = m%\ 



i™3/2 



nM n ~ 3 



\<t>\ 



2(n-l) 



(53) 



]\,{2(n-3) 



where is the mass of the AD field which is of order of 
the gravitino mass 771,3/2 in the gravity- mediation scenario 
and a m is a complex-valued coefficient of order unity. 
Here we have also included the Hubble-induced mass 
term with real- valued O(l) coefficient ch [13 , which 
is assumed to be positive here. In addition, the finite- 
temperature correction is significant as noted in the pre- 
vious section. This is represented as 



Kh(0)= £ c k f k \<j>\ 2 T 2 

fk\4>\<T 



( ^ 



(54) 

where and a is real and positive coefficient of order 
unity and fk is an appropriate coupling constant. Thus 
the potential of the AD field is given by 



V{</>) = Vs(4>) + Va l {<f>). 



(55) 



In the early Universe, the Hubble parameter is large and 
the AD field sits at the temporal minimum determined 
by the balance between the Hubble mass term and the 
nonrcnormalizable term, 



HM r 



_ 3 \l/(n-2) 



(56) 



The value of \<f>\ changes adiabatically as H decreases. 
The AD field begins the coherent oscillation when the 
Hubble parameter falls below the value H osc defined as 



rr2 

osc 



E 

fk\4>\<T 



Ckf k T 2 



a «s(T) 2 -J. (57) 



Because all contents of the Universe are diluted by 
the late-time entropy production, we need large baryon 
asymmetry enough to survive the dilution. Such a large 
initial baryon asymmetry can be generated by consider- 
ing the Affleck-Dine (AD) mechanism [46|, |47j- Actually 
the AD mechanism may work well even if the late-time 
entropy production dilutes the baryon number [48j . The 
MSSM contains many scalar fields which are scalar su- 
perpartners of the SM fermions, and there exist flat direc- 
tions in the scalar potential which do not feel the poten- 
tial in the SUSY limit at the renormalizablc level. Such 
a flat direction is called the AD field and is parametrized 
by <f>. The potential of the AD field is lifted up by the 
nonrcnormalizable superpotential 



W-nr 



nM 



n— 3 



with n > 4, 



In our scenario, it is assumed that the oscillation of 
the AD field is caused by the thermal-log potential, i.e. 
H 2 SC ~ ttgT 4 /|0| 2 , and the oscillation starts before the 
decay of the inflaton, _ff osc > T/, where T] is the decay 
rate of the inflaton. In this period, the temperature of 
the Universe is given by T w (T^HMp) 1 / 4 . Since the re- 
sultant baryon asymmetry is extremely small in the case 
of n = 4 [48|], we consider the case of n = 6. The baryon 
number density at the time when the AD field starts to 
oscillate is ns(t sc) ~ Scpm 3 / 2 \4'osc\ 2 , where <5cp is an 
effective CP phase which is assumed to be 0(1). After 
the reheating caused by the inflaton decay, the saxion 
starts to oscillate at t = ij. Since the saxion produces 
the present entropy in our model, the baryon-to-entropy 
ratio is fixed at that time. Using \(f> OS c\ — {HoscM 3 ) 1 / 4 
and H osc ~ (a 2 T|AfpA/ _3/2 ) 2/3 , the baryon-to-entropy 



(52) ratio after the decay of the saxion is estimated as 
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n B _ST a n B {U) _ 3T a n B (t osc ) ( a(t osc )\ 3 _ iT a n B (t osc ) ( T/ \ 2 / m 3 / 2 x 



s 4 Pa {ti) 4 Pa {ti) \ a(ti) ) 4 p a (ti) \H osc J \ Tj 

To \ ( m 3/2 \V 2 fl0 11 GcV\ fa s M P \ 2 ( M 

100A/p 



q in -i c ( T a \ / m 3/2 x 1/2 /10"GeV\ 

i3xl ° 5cp (iMcV ) [vEv) {—r^J 



Hence, the desired baryon asymmetry, n B /s ~ 10~ 10 , is 
generated by choosing the cutoff scale M appropriately. 
Notice that the baryonic isocurvature perturbation is suf- 
ficiently small and satisfies the observational bound jij. 

Let us check the consistency of the AD mechanism 
with present inflation model. During inflation, the AD 
field has VEV of ~ (HrM 3 ) 1 / 4 , where H I = nf 2 a /^M P 
is the Hubble scale during inflation. Thus the superpo- 
tcntial (|52[) induces an effective constant term given by 
W c ff ~ (iJ/M) 3 / 2 /6. Numerically, this is comparable to 
the constant term Wo = m 3 / 2 Mp for the above param- 
eter choices. Thus the inflaton dynamics is not affected 
by the presence of the large amplitude of the AD field. 

VIII. A VARIANT MODEL : SMOOTH-HYBRID 
INFLATION FROM A SUSY AXION MODEL 

So far we have analyzed the SUSY axion model based 
on the superpotential ([T]) and shown that it causes hybrid 
inflation. In this section we consider a variant type of the 
SUSY axion model. 

Let us take the superpotential 

W = S (iSS) - M 2 ) + A™, (58) 

where M is a cutoff scale and n > 2. In addition to the 
PQ symmetry and i?-symmetry whose charges are given 
in Table HI this superpotential also has a discrete symme- 
try Z n under which * has a charge +1 and others have 
zero. This has a flat direction along = (^M™" 1 ) 2 /' 1 , 
and it is stabilized by the SUSY-breaking masses similar 
to the model of (p}. Thus the minimum of the potential 
lies at |*| - |*| - f a = (//A/™- 1 ) 1 /™, where the PQ 
symmetry is spontaneously broken with scale of f a . 

The superpotential (|55|) . which provides one of the 
SUSY axion models, coincides with that causing a so- 
called smooth- hybrid inflation [4^, [5(J. Although the 
post-inflationary saxion dynamics is similar to the previ- 
ous model, there is a significant difference between hybrid 
and smooth-hybrid inflation. In the smooth-hybrid infla- 
tion model, the *(*) has a nonzero VEV during infla- 
tion, hence the PQ symmetry is already broken. There- 
fore, no topological defects are formed after inflation and 
we do not need to worry about the possibly harmful do- 
main wall problem. The scalar spectral index is pre- 
dicted to be around 0.97 without introducing nonminimal 
Kahlcr potentials. 



Instead, since the PQ symmetry is broken and the ax- 
ion obtains a quantum fluctuation during inflation, it 
may generate a large cold dark matter (CDM) isocur- 
vature perturbation. The amplitude of the CDM isocur- 
vature perturbation is given by 

Sm ^^^-, (59) 

where f2 m is the DM density parameter and ipN denotes 
the field value of the PQ scalar when the cosmological 
scales exit the horizon. The axion density parameter is 
evaluated in a way similar to (|47p , but in the present case 
the initial misalignment angle 9i can be chosen arbitrarily 
since the PQ symmetry is already broken during inflation 
and a region with angle Bi is expanded to cover the whole 
observable region of the Universe. Thus we have 

a,^o.o3r^(^)( T ^)V (60) 

By demanding that the correct magnitude of the density 
perturbation and the scalar spectral index of ~ 0.97 are 
reproduced, we find fi~3x 10 14 GeV and M - 10 15 GeV 
for n = 2. The situation does not change much for n > 2. 
Thus we obtain the PQ scale as / ~5x 10 14 GeV. Nu- 
merically, tpN is slightly smaller than f a . The inflation 
scale is calculated as Hi = p 2 /(V3M P ) - 10 10 GeV. Sub- 
stituting these parameters, we obtain S m ~ 10~ 6 for 
8i ~ 0.1 which is close to the observational bound on 
the isocurvature perturbation Q. In this case the axion 
cannot be the dominant component of DM. 

IX. CONCLUSIONS 

We have shown that an inflation naturally takes place 
in the framework a SUSY axion model. Identifying the 
PQ scalar fields as the waterfall fields, the hybrid in- 
flation is realized. In this case, the PQ phase transi- 
tion takes place at the end of inflation. Note that the 
isocurvature perturbation does not arise since the PQ 
symmetry is restored during inflation. In order for the 
inflation to account for the observed density perturba- 
tion, the PQ symmetry-breaking scale must be around 
10 15 GeV. In addition, the observed value of the spectral 
index can be reproduced by introducing a nonminimal 
Kahlcr potential. Considering the post-inflationary dy- 
namics, we have found that the saxion begins to oscil- 
late with large initial amplitude and its decay produces a 
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huge entropy. Thanks to this late-time entropy produc- 
tion process, the axion is diluted and its present density 
becomes consistent with that of DM. Simultaneously, the 
gravitinos [HI, HH and axinos [4(| are also diluted by the 
late-time entropy production, and hence the cosmolog- 
ical problems of gravitino or axino do not arise in the 
present scenario [a, [5l|. On the other hand, the pre- 
existing baryon asymmetry is also diluted, so we need 
an initial baryon asymmetry large enough to survive the 
dilution. We have pointed out that the desired baryon 
asymmetry is obtained through the AD mechanism. Fi- 
nally, topological defects formed at the PQ phase tran- 
sition or the QCD phase transition disappear without 
causing any cosmological problems if the color anomaly 
number is equal to 1. Therefore, the SUSY axion model 
studied in this paper not only provides a solution to the 
strong CP and the hierarchy problems but also accounts 
for inflation and the DM of the Universe. 

We have also shown that a slight modification of 
the SUSY axion model results in a so-called smooth- 



hybrid inflation. In this case no topological defects are 
formed and the scalar spectral index is consistent with 
the WMAP observation naturally. Instead, it predicts a 
large isocurvature perturbation close to the present ob- 
servational bound. Or it may be discovered by the future 
cosmological measurements. 



Acknowledgment 

K.N. would like to thank Fuminobu Takahashi for 
useful discussion. This work is supported by Grant- 
in-Aid for Scientific research from the Ministry of Ed- 
ucation, Science, Sports, and Culture (MEXT), Japan, 
No. 14102004 (M.K.), No. 21111006 (M.K. and K.N.), 
No. 22244030 (K.N.) and also by World Premier In- 
ternational Research Center Initiative (WPI Initiative), 
MEXT, Japan. 



[1] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 
(1977). 

[2] For reviews, see J. E. Kim, Phys. Rept. 150, 1 (1987); 

J. E. Kim and G. Carosi, Rev. Mod. Phys. 82, 557 (2 010). 
[3] For a review, see S. P. Martin, arXiv:hep-ph/9709356 
[4] E. Komatsu et al. [ WMAP Collaboration J, Astrophys. J. 

Suppl. 192, 18 (2011). [arXiv: 1001. 45381 [astro-ph.CO]]. 
[5] M. Kawasaki, N. Kitajima and K. Nakayama, Phys. Rev. 

D 82, 123531 (2010) [arXiv:1008.50l3l [hep-ph]]. 
[6] P. J. Steinhardt and M. S. Turner, Phys. Lett. B 129, 51 

(1983) ; G. Lazarides, C. Panagiotakopoulos and Q. Shan, 
Phys. Lett. B 192, 323 (1987); G. Lazarides, R. K. Schae- 
fer, D. Seckel and Q. Shan, Nucl. Phys. B 346, 193 
(1990); 

[7] M. Kawasaki, T. Moroi and T. Yanagida, Phys. Lett. B 

383, 313 (1996). 
[8] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shif- 

man, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 

166, 493 (1980). 
[9] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B 

104, 199 (1981); A. R. Zhitnitsky, Sov. J. Nucl. Phys. 

31, 260 (1980) [Yad. Fiz. 31, 497 (1980)]. 
[10] T. Kugo, I. Ojima, T. Yanagida, Phys. Lett. B135, 402 

(1984) . 

[11] T. As aka, M. Yamagu chi, Phys. Lett. B437, 51-61 
(1998) [hep-ph/98054491 ; Phys. Rev. D 59, 125003 (1999) 
arXiv: hep-ph /98 1 1451J . 

[12] N. Abe , T. Moroi, M Yamaguchi, JHEP 0201, 010 

(2002) . jhep-ph/0111155| . 

[13] T. Bank s, M. Dine, M. Gr aesser, Phys. Rev. D68, 075011 

(2003) . hep-ph/0210256| . 

[14] S. Nakamura, K. -i. Okumura, M. Yamaguchi, Phys. Rev. 

D77, 115027 (2008). i arXiv:0803.3"725l [hep-ph]] . 
[15] L. M. Carpenter, M. Dine, G. Festuccia, L. Ubaldi, Phys. 

Rev. D80, 125023 (2009). arXiv:0906.5015 [hep-th]]. 
[16] K. S. Jeong, M. Yamaguchi, arXiv: 1102.3301 [hep-ph]]. 
[17] T. Higaki and R. Kitano. larXiv: 1104.01701 [hep-ph], 
[18] E. J. Chun and A. Lukas, Phys. Lett. B 357, 43 (1995). 



[19] M. Kawasaki, K. Nakayama and M. Senami, JCAP 0803, 
009 (2008). 

[20] E. J. Copeland, A. R. Liddle, D. H. Lyth, E. D. Stew- 
art and D. Wands, Phys. Rev. D 49, 6410 (1994); 
G. R. Dvali, Q. Shan and R. K. Schaefer, Phys. Rev. 
Lett. 73, 1886 (1994); G. Lazarides, R. K. Schaefer and 
Q. Shan, Phys. Rev. D 56, 1324 (1997); A. D. Linde and 
A. Riotto, Phys. Rev. D 56, 1841 (1997); G. R. Dvali, 
G. Lazarides and Q. Shan, Phys. Lett. B 424, 259 
(1998); V. N. Senoguz and Q. Shan, Phys. Rev. D 71, 
043514 (2005); R. Jeannerot and M. Postma, JHEP 
0505, 071 (2005); M. Bastero-Gil, S. F. King and 
Q. Shan, Phys. Lett. B 651, 345 (2007); M. ur Rehman, 
V. N. Senoguz and Q. Shan, Phys. Rev. D 75, 043522 
(2007); M. U. Rehman, Q. Shan and J. R. Wickman, 
Phys. Lett. B 683, 191 (2010); C. Pallis, JCAP 0904, 
024 (2009). 

[21] K. Nakayama, F. Takahashi, T. T. Yanagida, JCAP 

1012, 010 (2010). 
[22] M. Bastero-Gil and S. F. King, Phys. Lett. B 423, 27 

(1998). 

[23] G. Lazarides and C. Pallis. larXiv:1007.1558l [hep-ph], 
[24] S. R. Coleman and E. J. Weinberg, Phys. Rev. D 7, 1888 

(1973). 

[25] A. R. Liddle and D. H. Lyth, "Cosmological Inflation 
and Large-scale Structure, " (Cambridge University Press, 
Cambridge, England, 2000). 

[26] A. D. Linde, Phys. Rev. D53, 4129-4132 (1996) 
[hep-th/9601083| . 

[27] A. Anisimov and M. Dine. Nucl. Phys. B 619, 729 (2001). 

[28] Y. I. Izotov, T. X. Thuan, Astrophys. J. 710, L67- 
L71 (2010). [arXiv: 1001. 44401 [astro-ph.CO]]; E. Aver, 
K. A. Olive, E. D. Skillman, JCAP 1005, 003 (2010). 
[arXiv:1001. 52181 [astro-ph.CO]]; E. Aver, K. A. Olive, 
E. D. Skillman, [arXiv:1012.2385l [astro- ph.CO]]. 

[29] J. Dunkley, R. Hlozek, J. Sievers et al., |arXiv:1009.'0866l 
[astro-ph.CO]]. 

[30] M. Kawasaki, K. Kohri and N. Sugiyama, Phys. Rev. 



12 



Lett. 82, 4168 (1999); Phys. Rev. D 62, 023506 

(2000); S. Hannestad, Phys. Rev. D 70, 043506 (2004); 

K. Ichikawa, M. Kawasaki and F. Takahashi, Phys. Rev. 

D 72, 043522 (2005). 

[31] J. Hasenkamp and J. Kersten, arXiv:1103.6193 [hep-ph]. 
[32] D. J. Gross, R. D. Pisarski and L. G. Yaffe, Rev. Mod. 

Phys. 53, 43 (1981). 
[33] M. S. Turner, Phys. Rev. D 33, 889 (1986). 
[34] M. Bolz, A. Brandenburg and W. Buchmuller, Nucl. 

Phys. B 606, 518 (2001) [Erratum-ibid. B 790, 336 

(2008)]; J. Pradler and F. D. Steffen, Phys. Rev. D 75, 

023509 (2007); Phys. Lett. B 648, 224 (2007). 
[35] M. Kawasaki, F. Takahashi and T. T. Yanagida, Phys. 

Lett. B 638, 8 (2006); Phys. Rev. D 74, 043519 (2006); 

T. Asaka, S. Nakamura and M. Yamaguchi, Phys. Rev. 

D 74, 023520 (2006); M. Endo, K. Hamaguchi and 

F. Takahashi, Phys. Rev. D 74, 023531 (2006); M. Endo, 

F. Takahashi and T. T. Yanagida, Phys. Lett. B 658, 236 

(2008); Phys. Rev. D 76, 083509 (2007). 
[36] M. Kawasaki, K. Kohri and T. Moroi, Phys. Lett. B 625, 

7 (2005); Phys. Rev. D 71, 083502 (2005); M. Kawasaki, 

K. Kohri, T. Moroi and A. Yotsuyanagi, Phys. Rev. D 

78, 065011 (2008). 
[37] K. Jedamzik, Phys. Rev. D 74, 103509 (2006). 
[38] K. Rajagopal, M. S. Turner and F. Wilczek, Nucl. Phys. 

B 358, 447 (1991). 
[39] T. Goto and M. Yamaguchi, Phys. Lett. B 276, 103 

(1992). 



[40] L. Covi, H. B. Kim, J. E. Kim and L. Roszkowski, JHEP 
0105, 033 (2001); A. Brandenburg and F. D. Steffen, 
JCAP 0408, 008 (2004) ; A. Strumia, JHEP 1006, 036 
(2010) |arXiv: 1003.58471 [hep-ph]]. 

[41] C. Cheung, G. Elor and L. J. Hall. larXiv:1104.0692l [hep- 
ph]. 

[42] A. Vilenkin and A. E. Everett, Phys. Rev. Lett. 48, 1867 
(1982). 

[43] M. Yamaguchi, M. Kawasaki and J. Yokoyama, Phys. 
Rev. Lett. 82, 4578 (1999). 

[44] T. Hiramatsu, M. Kawasaki, T. Sekiguchi, M. Yam- 
aguchi, J. 'i. Yokoyama, (arXiv: 1 012. 5502 [hep-ph]]. 

[45] H. Georgi and M. B. Wise, Phys. Lett. B 116, 123 (1982); 
J. E. Kim and H. S. Song, Phys. Rev. Lett. 55, 677 
(1985). 

[46] I. Affleck and M. Dine, Nucl. Phys. B 249, 361 (1985). 
[47] M. Dine, L. Randall and S. D. Thomas, Nucl. Phys. B 

458, 291 (1996). 
[48] M. Kawasaki and K. Nakayama, Phys. Rev. D 76, 043502 

(2007) . 

[49] G. Lazarides, C. Panagiotakopoulos, Phys. Rev. D52, 

R559-563 (1995) |hep-ph/9506"325] . 
[50] M. Ya maguchi, J. 'i. Yok oyama, Phys. Rev. D70, 023513 

(2004) |hep-ph/0402282] . 
[51] M. Kawasaki and K. Nakayama, Phys. Rev. D 77, 123524 

(2008) . 



